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1.  Introduction 

Many  of  the  most  important  and  influential  recent  works  which 
dynamic  meteorologists  have  accepted  as  belonging  in  their  field,  and 
which  have  been  published  in  meteorological  journals,  have  dealt  with 
the  behavior  of  a  gas  of  uniform  composition.  By  contrast,  the  most 
visible  features  of  the  real  atmosphere,  as  is  especially  evident  when 
one  looks  at  satellite  photographs  of  the  earth,  are  the  clouds,  whose 
variable  occurrence  in  space  and  time  demands  a  fluid  of  varying 
composition.  The  clouds  are  not  mere  atmospheric  ornaments;  the  fact 
that  they  are  visible  implies  that  they  reflect  much  of  the  sun's 
radiation  which  would  otherwise  penetrate  to  lower  levels  in  the 
atmosphere  and  would  heat  the  ocean  and  land  surfaces  underneath. 
Moreover,  the  clouds,  and  the  water  vapor  which  must  be  present  if  the 
clouds  are  not  to  evaporate  immediately,  are  strong  absorbers  and 
emitters  of  long-wave  radiation,  while  the  gain  and  release  of  latent 
heat,  which  must  accompany  changes  in  pressure  if  the  air  is  to  remain 
saturated,  give  moist  air  a  different  thermodynamical  behavior  from  dry 
air.  It  is  therefore  relevant  to  pose  the  question,  "To  what  extent  is 
the  dry  atmosphere  which  appears  in  so  many  theoretical  studies  an 
acceptable  model  of  an  atmosphere  where  water  in  its  gaseous,  liquid, 
and  solid  phases  occurs  in  varying  concentrations?" 

One  can  of  course  get  some  idea  of  the  answer  to  this  question  by 
comparing  the  results  of  these  theoretical  studies  with  observations  of 
the  real  atmosphere,  or  with  other  theoretical  studies  where  the 
presence  of  water  enters  explicitly.  The  latter  studies  logically 
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inelude  the  recent  work  with  large  operational  global  circulation 
models.  However,  dynamlcists  often  attack,  questions  of  this  sort  by 
constructing  less  complicated  models,  in  which  some  real  properties  and 

processes  deemed  to  have  no  bearing  on  these  questions  are  purposely 

omitted.  The  most  extreme  models  of  this  sort  are  the  low-order  models, 
formulated  with  a  minimal  number  of  dependent  variables,  where  all  but 
the  essential  features  have  been  eliminated,  and  where  care  must  often 
be  taken  not  to  eliminate  some  of  the  essentials  as  well. 

Most  low-order  models  of  circulations  of  global  scale  have  been  dry 
models.  It  is  easy  to  see  why  this  ha6  been  so.  In  constructing  a 

low-order  model  one  usually  strives  for  simplicity.  Aside  from  the 
obvious  fact  that  there  are  more  basic  dependent  variables  to  handle  in 
a  moist  model,  both  the  radiative  and  the  thermodynamic  processes  which 
characterize  a  moist  atmosphere  introduce  complicated  nonlinear  terns 
Into  the  governing  dynamical  equations.  Low-order  models  are  most 
frequently  formulated  in  terms  of  spatially  orthogonal  functions,  since, 
when  sparsely  distributed  grid  points  are  used,  horizontal 

finite  differences  do  not  generally  afford  good  approximations  to 
horizontal  derivatives.  Simple  products  are  easy  to  evaluate  in 
orthogonal-function  s pace,  but  more  complicated  nonlinear  functions  are 
no  t . 

In  large  moist  models  which  use  orthogonal  functions,  this 
difficulty  is  generally  circumvented  by  transforming  the  variables  from 
orthogonal-function  space  to  grid-point  space  at  each  time  step, 
evaluating  the  cumbersome  nonlinear  terms  at  each  grid  point,  and  then 
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transforming  back  to  orthogonal-function  space,  where  the  horizontal 
derivatives  are  evaluated.  This  procedure  Is  also  possible  in  a 
low-order  moist  model,  hut  It  introduces  other  dlf f iculltes.  Because 
a  global  field  of  water  vapor  content  possesses  horlzo  ital  variations 
comparable  to  its  mean  value,  the  representation  of  a  realistic  field  in 
terms  of  a  small  number  of  orthogonal  function,  with  coefficients  whi^h 
are  optimal  for  the  bulk  of  the  atmosphere,  is  likely  to  vield 
supersaturation  when  evaluated  at  some  tropical  grid  points,  and 
negative  moisture  content  at  some  polar  grid  points.  Neither  outcome  is 
tolerable. 

In  a  Final  Report  entitled  "The  feasibility  of  a  low-order  model  of 
a  moist  general  circulation",  prepared  under  a  previous  contract  with 
the  Air  Force  Geophysics  Laboratory,  we  noted  these  difficulties,  and 
described  possible  methods  for  overcoming  them.  In  brief,  we  suggested 
eliminating  the  possibility  of  supersauration  by  using  a  measure  of 
total  water  content  rather  than  water  vapor  content  as  a  basic  variable, 
and  introducing  an  auxiliary  equation  to  define  the  water  vapor  content 
in  terms  of  the  actual  total  water  content  and  the  water  vapor  content 
needed  to  produce  saturation.  We  proposed  to  avoid  the  likelihood  of 
negative  moisture  content  by  using  the  total  dew  point,  i.e.,  the  dew 
point  which  would  result  If  all  liquid  water  were  converted  to  vapor, 
rather  than  the  total-water  mixing  ratio  as  a  basic  variable;  the 


variations  of  total  dew  point,  like  those  of  temperature,  are  generally 
small  compared  to  mean  values. 
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We  then  proceeded  to  describe  in  detail  a  low-order  model  of  a 
moist  general  circulation  incorporating  the  suggested  features.  We 
described  some  preliminary  tests  of  the  model  for  the  special  case  of 
horizontally  uniform  solar  heating.  Tests  with  variable  solar  heating 
were  yet  to  be  performed. 

During  the  two  years  since  the  appearance  of  that  report  we  have 
made  some  substantial  changes  in  the  model.  We  have  done  so  for  a 
variety  of  reasons.  In  some  instances  the  possibility  of  making  the 
changes  simply  had  not  occurred  to  us.  In  other  instances  we  had  wished 
all  along  to  make  the  changes,  but  had  previously  been  unable  to  find  a 
satisfactory  way  to  do  so.  In  still  other  instances  we  had  not 
considered  It  desirable  or  necessary  to  make  the  changes,  and  were  led 
to  do  so  only  after  discovering  that  the  model  as  originally  formulated 
would  not  perform  properly.  Examples  of  changes  made  for  these  various 
reasons  are,  respectively,  the  use  of  vertical  averging  to  simplify  the 
vertical  structure  of  the  model,  the  inclusion  of  the  advection  of 
temperature  and  moisture  by  the  divergent  part  of  the  wind  field,  and 
the  use  of  a  grid  of  nine  instead  of  sixteen  points  for  evaluating  the 
nonlinear  radiative  and  thermodynamic  terms. 

Throughout  the  development  of  the  model  our  goal  has  been 
simplicity.  By  this  we  mean  simplicity  in  concept,  simplicity  in 
formulation,  and  simplicity  in  execution.  It  is  not  certain  that  we 
have  achieved  this  goal;  our  model  is  much  more  complicated  than  a  dry 
model  with  similar  vertical  and  horizontal  resolution. 


Because  of  our  attempts  to  simplify,  we  cannot  expect  the  model  to 
duplicate  the  behavior  of  the  earth's  atmosphere  in  more  than  a 
qualitative  sense.  We  might  better  regard  our  model  as  a  model  of  some 

hypothetical  planet's  atmosphere,  possessing  one  constituent  which 

occurs  in  two  phases.  However,  when  numerical  values  are  required,  we 
have  chosen  them  to  be  representative  of  the  earth's  atmosphere,  with 
water  as  the  varible  constituent. 

Because  we  anticipate  performing  comparative  numerical  studies  in 
which  some  constants  or  initial  values  are  altered  by  very  small 

amounts,  we  have  made  a  point  of  eliminating  all  discontinuities  in  the 
formulation  of  the  equations.  In  particular,  we  have  used  an 

unconventional  formula  for  the  water-vapor  mixing  ratio  in  terms  of 
temperature  and  total  water  content. 

In  the  following  sections  we  shall  present  a  detailed  description 
of  the  construction  of  the  model,  including  those  features  which  we  also 
described  in  the  previous  final  report.  We  shall  then  examine  some 
preliminary  numerical  solutions.  Finally,  we  shall  consider  briefly 
some  further  modifications  of  the  model  which  should  perhaps  eventually 
be  made. 


2.  The  continuous  equations 

Our  model  atmosphere  will  be  composed  of  a  mixture  of  dry  air, 
water  vapor,  and  liquid  water  in  variable  proportions.  It  will  be  of 


infinite  upward  and  west-east  extent,  and  will  be  bounded  on  the  south 
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and  north  by  frictionless  vertical  walls.  Under  the  atmosphere  will  be 
an  ocean  of  uniform  composition  and  limited  vertical  extent.  The 
atmosphere  will  be  governed  by  baroclinic,  geostrophic,  beta-plane 
dynamics,  while  the  ocean  will  exchange  heat  and  water  with  the 
atmosphere.  The  sole  external  driving  force  will  be  incoming  solar 
radiation;  dissipation  will  be  both  thermal  and  mechanical.  We  shall 
not  consider  the  possible  presence  of  ice,  even  at  subfreezing 
temperatures,  nor  shall  we  allow  ocean  currents  to  develop. 

We  shall  begin  b>  formulating  the  equations  for  a  vertically  and 
horizontally  continuous  atmosphere-ocean  system.  We  shall  then  reduce 
the  atmospheric  equations  to  those  of  a  two-layer  model,  and  the  oceanic 
equation  to  that  of  a  single-layer  ocean,  by  vertical  averaging. 
Finally  wc  shall  introduce  a  truncated  set  of  orthogonal  functions,  in 
terms  of  which  each  horizontally  continuous  field  will  be  expressed. 

Our  twsic  equations  wilL  contain  the  constants 

width  of  atmospheric  channel,  divided  by  7?  , 
sea-Ievei  pressure, 

pressure  at  bottom  of  oceanic  layer, 

Coriolis  parameter  in  middle  of  channel, 
northward  derivative  of  Coriolis  parameter, 
specific  heat  of  air  at  constant  pressure, 
gas  constant  for  air, 
latent  heat  of  condensation, 
specific  heat  of  water, 
the  Independent  variables 
t  time. 
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x 

y 

P 

and  the 

V 

X 

'-V 

CO 

"t 

w 

/S/ 

s 


eastward  distance, 
northward  distance, 
pressure , 

basic  dependent  variables 
stream  function, 

velocity  potential  (with  sign  changed), 
individual  pressure  change, 
atmospheric  temperature, 
total  dew  point, 
oceanic  temperature. 


The  horizontal  wind  field  will  be  the  sum  of  a  rotational  nondlvergont 
part,  derivable  from  Vf'  ,  and  an  irrotational  divergent  part,  derivable 
from  X  •  To  this  list  we  shall  add  the  auxiliary  variables 

u  saturation  mixing  ratio  at  temperature  T  and  pressure  p  , 
v  water  vapor  mixing  ratio. 


w  total-water  mixing  ratio, 

s  saturation  mixing  ratio  at  temperature  S  and  pressure  p 


We  have  placed  a  tilde  ( 'r "')  over  each  dependent  variable  because  we  shall 
be  using  the  same  symbols,  without  tildes,  for  the  values  of  these 
variables  at  the  atmosphere-ocean  interface. 

The  basic  prognostic  equations  will  be  the  vortlcity  equation 


—  J($,  Vay)  -f 


(O 


-8- 


the  atmospheric  thermodynamic  equation 

A  ( Cp  T+L.M'^y4'*-~RT(-oyp»  +  H  J  (2) 

the  water-content  equation 

=  &  ,  (3) 

and  the  oceanic  thermodynamic  equation 

S)/Jt  =  E  .  <-'■> 

Here  .1  denotes  a  Jacobian  with  respect  to  x  and  y  ,  while  the 
diabatic  terms  F,  H,  C,  and  E  denote  respectively  the  curl  of  the 
viscous  drag,  the  diabatic  atmospheric  heating  per  unit  mass  (including 
chc  effect  of  evaporation  from  the  ocean  surface,  but  not  conder  '-at  loti 
within  the  atmosphere,  which  increases  T  but  decreases  v  ),  the  gain 
or  loss  of  water  by  evaporation  or  precipitation,  and  the  oceanic 
heating  per  unit  mass.  In  keeping  with  the  geostrophic  simplification 
we  have  omitted  the  advection  of  vorticity  by  the  divergent  part  of  the 
wind,  bv;t  the  complete  advection  of  enthalpy  and  water  are  implicitly 
included  in  the  time  derivatives  in  (2)  and  (3).  These  equations  are 
to  be  accompanied  by  the  diagnostic  continuity  and  thermal-wind 
equations 

3w/ap  =  -V  DC 


(5) 


Eq.  (6)  implies  that  Che  horizontal  average  of  *  increases  with 

— 

elevation,  and  that,  aside  from  a  constant  factor,  1^7  is  the  i.s&baric 
height.  For  boundary  conditions  we  shall  let  LaJ  =  0  when  p  =  o Q  or 

0  . 

The  auxiliary  dependent  variables  are  related  to  the  basic 
variables  by  the  formulas 


Here  c*  is  to  be  chosen  to  produce  the  proper  saturation  mixing  ratio 
at  some  standard  temperature  and  pressure,  and  the  choice  of  j  (and 
indeed  the  choice  of  Eq.  (10)  to  define  T  )  is  somewhat  arbitrary. 
Eqs.  ( 7 )  — ( 9)  may  be  derived  from  the  C lausius-C lapeyron  equation  if  the 
factor  T  in  that  equation  is  first  replaced  by  the  product  T*T  , 
where  T*  is  a  standard  temperature  typical  of  the  atmosphere,  say  273 
K.  This  approximation  makes  (8)— ( 10)  algebraic  instead  of 
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transoendental.  An  appropriate  value  for  the  exponent  r  is  about 
20.1',  and  by  letting  JtA  *20  exactly  we  can  effect  a  considerable  saving 

/V 

in  computation.  Eq.  (10)  makes  the  liquid  water  content  //v-  - /sT  small 

Aw 

when  the  degree  of  subsaturation  /A  -  rjr  is  large,  and  vice  versa. 
Choosing  Y  -  1/4  makes  the  relative  humidity  at  j  /CC  equal  to  80  percent 
when  T  and  W  are  equal;  the  remaining  20  per  cent  of  the  water  is  in 
the  form  of  clouds.  Choosing  V  =0  would  reduce  (10)  to  the 
conventional  assumption  that  there  is  no  supersaturation,  and  no  clouds 
with  subsaturation,  but  it  would  also  introduce  the  usual  discontinuity 
at  the  saturation  point;  we  therefore  prefer  a  posit  i  ve  y  .  Eqs. 
(l)-(10)  form  a  closed  system  in  the  ten  basic  and  auxiliary  dependent 
variables. 

Our  low-order  model  is  to  consist  of  a  simpler  system  of  equations, 
obtained  by  introducing  simplifying  assumptions  into  (1)— (10).  However, 
before  we  alter  the  meaning  of  (i)-(10)  at  alL,  we  can  realize  some 

simpiif iction  by  expressing  the  constants  and  variables  in  dimensionless 

X  X, 

form.  We  can  effectively  do  this  by  choosing  1/f,  D,  pfi  ,  and  D  f  /R 
as  the  units  in  which  time,  distance,  pressure,  and  temperature  are 
measured.  The  numerical  values  of  f,  D,  p  ,  and  K  are  then  unity, 
while  the  numerical  values  of  the  remaining  constants  are  the  ratios  of 
these  constants  to  the  combinations  of  f,  1),  p  ,  and  R  with  the 
appropriate  dimensions.  In  the  equations  which  follow  we  shall  assume, 
unless  we  state  otherwise,  that  all  quantities  have  been  made 
dimensionless,  and  we  shall  omit  any  factors  whose  numerical  values  must 


be  unity 


3.  Vertical  simplification 


in  many  two-  level  or  two-layer  geostrophic  models  the  stream 
function  is  specified  at  each  of  two  levels,  or  averaged  through  each  of 
two  layers,  while  the  temperature  is  specified  for  only  one  level  or 
layer.  The  temperature  is  then  identified  through  the  thermal -  wind 
equation  with  the  difference  of  the  stream  functions.  We  shall  use  the 
two-stream-function  one-temperature  format  in  our  model.  It  is 
consistent  with  this  formulation  to  specify  the  total  dew  point  for  only 
one  level  or  layer. 

Our  procedure  for  reducing  Eqs.  (1)-(10)  to  a  two-layer  model  will 
be  dictated  by  a  consideration  of  the  source  and  sink  terms  F,  H, 
and  G  .  It  is  not  at  all  obvious  how  these  should  be  formulated  at 
individual  levels  within  the  atmosphere.  However,  if  horizontal 
dissipative  exchanges  of  momentum,  energy,  and  water  in  the  atmosphere 
are  considered  negligible,  the  vertical  averages  of  F,  H,  and  G 
simply  represent  net  exchanges  between  the  atmosphere  and  the  ocean, 
and.  In  the  case  of  energy,  between  the  atmosphere  and  outer  space. 
Accordingly,  we  shall  construct  our  vertically  simplified  model  by 
averaging  Eqs.  (2)-(4)  through  the  depth  of  the  atmosphere  or  ocean. 
Eq.  (1)  must  be  averaged  separately  through  each  of  two  layers  In  order 
to  yield  two  stream  functions,  whence  the  exchange  of  momentum  between 
the  two  layers  will  also  enter  the  model. 

Because  u,  v,  and  w  normally  fall  :>ff  much  more  rapidly  with 
elevation  than  T  and  W  ,  we  cannot  choose  any  single  level  and  claim 
that  the  values  of  the  various  terms  in  (l)-(3)  at  this  level  are 
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representat 1 ve  of  vertical  averages.  We  shall  therefore  specify  the 

manner  in  which  each  dependent  variable  varies  with  p  ,  whereupon 
vertical  integration  will  be  feasible.  We  shall  find  it  convenient  to 
express  the  dependent  variables  in  terms  of  their  values  at  the 
atmosphere-ocean  interface;  these  will  be  denoted  by  symbols  without 
tildes,  and  will  become  our  new  dependent  variables. 

We  first  let  the  lapse  rate  of  atmospheric  temperature  with 
elevation  be  constant,  so  that 

— '  A 

T  '  T  p  ,  (in 

where  is  a  constant.  We  make  no  attempt  to  model  a  stratosphere. 

Since  =  0  would  imply  an  isothermal  stratification,  and  X  =  K  = 

^'cp  would  imply  a  dry-adiabat ic  lapse  rate,  X  should  be  between  0 
and  X  .  We  omit  ail  vertical  variations  of  oceanic  temperature, 
since  we  are  considering  only  a  shallow  layer,  so  that 

5  =  S. 


ft  follows  from  (6),  (7),  and  (9)  that 

Y  =  Y  +  T  (  I  -  p*)/x 

~  /*  x  *  i 

M  2  A*  Y 


(14) 
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S  =  s 


(13) 


Next  we  let  the  relative  humidity  v/u  be  constant  within  each 
vertical  column.  Then 


rr  -  /vr  f 


/*  >  - 


06) 


Eq.  (10),  being  homogeneous  in  u,  v,  and  w  , still  holds  when  the 
tildes  are  dropped,  so  that 


MT  -  _M/~  j? 


(17) 


whence,  from  (8), 


W  =  Wf> 


(18) 


In  order  that  remain  finite  as  the  top  the  atmosphere  is  approached, 

it  is  necessary  that  X  >  V  ,  where  V  *  '//*  =  0.05.  Finally  we  let 

A,  A# 

*  ,  like  'f'  ,  be  linear  in  p  ,so  that,  if  (5)  is  to  be  satisfied, 

rs* 

with  tu  =  0  when  p  =•  0  or  1, 


y.  =  dc  ( -  i  +  ( i  +  x)  p")/  y  j 


(19) 


(20) 
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We  may  treat  Eqs.  (ll)-(20)  as  approximations  to  be  used  only  tor 
vertical  averaging,  or  we  may  assume  that  K,  H,  C,  and  E  implicitly 
include  Internal  processes  which  6erve  to  maintain  the  prescribed 
vertical  distributions  against  the  disruptive  effects  of  advection  and 
other  processes,  just  as  vertical  motion  and  divergence  serve  to 
maintain  hydrostatic  and  geostrophic  equilibrium  against  the  disruptive 
effects  of  such  processes. 

We  shall  average  Eq.  (1)  separately  from  p  =  0  to  1  and  p  = 

t  i  f 

to  a  level  A  ;  this  is  equivalent  to  averaging  from  0  to  A  and  X 

I 

to  1  .  The  most  logical  value  of  X  would  probably  be  1/2,  but  it 
makes  the  equations  slightly  simpler  to  let  A  -  (  I  —  X)  ^  >  or,  f°r 

acceptable  values  of  A  >  about  1/3.  We  then  find  that  two  linear 
combinations  of  the  averaged  equations  are 


-V  X  +  -  F 


sv’t/h  =  •  Jit,  v't)  - t(t,vV) bn)j(T vt) 

-fiT/ s* +(n->.)vV -(.n)(F-F)  ,  ‘ 


where  the  .iouhle  and  single  bars  over  r*  denote  averages  from  0  to  1 
and  0  to  X  >  respectively. 

Averaging  Eq.  (2),  after  using  (10),  (7),  and  (3)  to  express  dv/dt 
in  terms  of  dT/dt  and  G  ,  we  obtain 
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^*7  t+,T))r  'b*  iy'v^)^r  7]C 
+((*.->)  -v(x-yW)rvV  +  (<+ x)(.+^(  a.  -  , 


where 


l_  (1*>)  (»+?•>) 


^  ^  aa 

}M  T 


(2-) 


and  the  partial  derivatives  J  v/Ju  and  3v/9w  are  to  be  evaluated 
from  (10).  Averaging  F.q.  (3),  and  then  using  (8)  to  express  w  in 
terms  of  W  ,  we  obtain 


x(l+^(^  -  -V  J(T,w)  »-M?w-vx 

-V(x-y)  W  V  X  +  A  I1'1'*’)  ~£r  G-  .  c 


Finally,  Eq.  (4)  becomes 

-  B/c  ,  (26) 

With  the  auxiliary  variables  u,  v,  w,  and  s  defined  by  (7)— (10), 
applied  at  p  =  1  ,  and  with  the  definition  (24),  Eqs.  (21  )— (23)  and 
( 2 5 )  —  ( 26 )  become  a  closed  sytem  in  the  basic  dependent  variables  y'  , 
~y C.  ,T,  W,  and  S 
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t.  Horizontal  simplification 

Low-order  models  of  dry  atmospheric  circulations  are  typified  by  .< 
two-layer  model  which  we  introduced  some  time  ago  (Lorenz,  196 'J)  to 
study  the  phenomenon  of  vacillation,  i.e.,  periodic  behavior  which 
cannot  be  reduced  to  steady  behavior  by  introducing  a  moving  coordinate 
system.  In  that  model  we  expanded  each  dependent  variable  in  a 
truncated  series  of  orthogonal  functions;  specifically,  in  a  somewhat 
different  notation,  we  let 
6 

A  -  £  Aj, 

«*o  >  v  ' 


where  A  stands  for  any  dependent  variable,  and 

' 

<£  ,  a  3.  y.  CoS  a  % 

f  ^  s  7. 

<-oS  ** 

r  2.  £>* 

$<,  ^  cos 

'I  lie  orthogonal  functions  were  chosen  so  that 

*  M 


(2K) 


^4.  --  -aK  5 


(29) 


where  *  0,  a  (  »  a  ^  »  5,  a  ^  •  1 ,  anci  a  6  ~  They 


satisfy  the  orthonorma i i ty  conditions 
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On) 


so  that,  if  C  is  any  quantity  to  be  expanded  in  a  truncated  series, 

,~n  tr 


c„--  •*"(  ( 

'A  ^  A 


(31) 


We  then  substituted  the  expansion  (27)  for  each  dependent  variable 
into  the  horizontally  continuous  equations.  Upon  expanding  the  left  and 
right  sides  of  the  resulting  equations  in  truncated  series,  using  (31), 

and  then  equating  the  coefficients  of  CD  , . . . ,  <£  ,  we  effectively 

O  x6 

converted  each  continuous  equation  ir.ro  seven  ordinary  differential 
equations;  these  became  the  equations  of  the  low-order  model. 

In  the  present  model  we  shall  express  the  dependent  variables  in 
terms  of  the  same  truncated  set  of  orthogonal  functions.  We  may,  in 

fact,  regard  the  present  model  as  an  extension  of  the  earlier  model  to  a 
moist  atmosphere. 

In  the  dry  model  the  right-hand  sides  of  the  continuous  equations 
contain  nonlinear  terms  of  the  form  J(A,B),  to  be  expanded  into 

orthogonal  functions.  According  to  (31),  if  C  =  -  .T(A,B), 


=  o 

Cl 

r  S 

bl\3]x. 

+  b  C  A;  8} 

Cpl 

--  5 

b  [a,s]3, 

t  ^  b  la,b7 

Cj 

-  5 

b  £A;8)ii 

■+  W  b 

Ch 

% 

b  L.a, 81^ 

+  Wb  Ca,8]sj 

cs 

=  % 

b  LA/  B^]  t  ( 

--  % 

t>  [A,B]  ,s 

<3  b  [A , v 

(32) 


J 


-18- 


whore  b  -  (16/2/15 TT ),  and  lA,B)mn  stands  for  AmBn  -  AnBm  .  The 
present  model  also  contains  terms  of  the  form  AB  .  From  (31),  if  C  = 
AB  , 

Co  r  r  Ak  B*  A 

V\  -  «  ) 

c,  --  (A, s)„  -  ((A,3),t-CA,e )„)/r* 

C,  --  (A,  8).»  -  ((A,b)„- (A,3)„V^ 

+  ((A,  8\,  +  (A,B)>  y  +  (a,  0),  ; 

C,  '  (A, 8)  OH  ~  '  (A,8).,)/l> 

C  s  =  (.A,  8)„,  -t(AyS)st  -Lw)»)in. 

C<,  r  (A,e\t  •  ((A,8)<t  +A,8,tAiax  +  A>B,yr» 

where  (A, R )nm  stands  for  AmBn  +  AnBm  .  Formulas  for  Cn  ,  where 
c  -  a-7b,  can  be  derived  from  (29)  and  (33)  by  noting  that 

7A-V6  «  l3-> 

As  already  noted,  to  evaluate  the  more  complicated  non 3 i near 
functions  appearing  in  the  source  and  sink  terms,  and  in  the 

mulst-thermodynamie  terms,  it  Is  preferable  to  transform  from  orthogonal 
functions  to  grid  points,  evaluate  the  nonlinear  terms  at  each  grid 

point,  and  then  transform  back.  We  have  chosen  a  set  of  nine  grid 

points,  at  the  intersections  of  the  lines  x  =  TT/b,  3’h'/b,  5  7T/t>  with 
the  lines  y  *  7f  lb,  3  7T /b,  5  / b.  Our  choice  of  nine  points  has 
followed  earlier  experimentation  with  a  grid  of  lb  points.  We 
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discovered,  for  the  case  where  there  were  no  variations  with  x  , 
that  unreasonably  high  temperatures  would  develop  if  the  number  of  grid 
points  In  the  y-directlon  exceeded  the  number  of  orthogonal  functions 
(here  Q,  3,  <£  (,)  depending  upon  y  alone.  This  difficulty 

disappeared  when  the  number  of  grid  points  was  reduced  to  the  number  of 
orthogonal  functions.  This  result  appears  contrary  to  our  experience 
with  transformations  between  orthogonal  functions  and  grid  points  when 
the  nonlinear  terms  are  quadratic,  and  apparently  results  from  the 
higher-degree  nonlinearity.  Thus  far  we  have  encountered  no 

difficulties  from  the  use  of  nine  rather  than  seven  grid  points,  and,  in 
any  event,  we  have  found  no  way  to  choose  a  reasonably  spaced  grid  of 
seven  points. 

Our  model  contains  no  prognostic  equation  for  X  .  To  integrate 
the  equation,  we  must  eliminate  D  T  /  5 1  from  (22)  and  (23),  and  solve 
the  resulting  diagnostic  equation  for  •  Although  some  of  the  terms 

must  be  evaluated  in  grid-point  space,  the  equation,  which  contains 
horizontal  derivatives  of  ^)C  ,  is  not  suitable  for  solution  in  grid- 

point  space,  and  special  care  must  be  taken  to  make  it  readily  solvable 
even  In  orthogonal-function  space. 

-a.  a. 

We  note  first  that  if  y  denotes  the  particular  inverse  of 
whose  horizontal  average  vanishes,  (22)  and  (23)  may  be  written  as 

?(T  -T.)/at  *  w'toc-x.)  +  V* 1 F  "  ,  (30 

x'  3T/>t  =  y’vT-  v  X  +2TvaX  +  h'  ,  OM 
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where  F'  and  H’  include  all  Che  terms  In  (22)  and  (23)  which  contain 
neither  nor  a  time  derivative.  It  is  easy  to  eliminate  2  T/3  t. 

Without  loss  of  generality  we  can  let  =  0.  According  to  (33)  and 

(34) ,  when  the  resulting  equation  is  transformed  to  orthogonal-function 

space,  terms  of  the  form  A")C  or  Ay  >  will  transform  to  linear 

expressions  in  the  unknowns  "]£f,  ,  whose  coefficients  are  linear 

functions  of  Afl,  ...,  A^.  However,  the  term  y-v  T.qyi  is  not  of  the 
proper  form.  Accordingly,  we  first  divide  (36)  by  Y'  ,  which 
fortunately  never  vanishes,  and  then  subtract  the  proper  multiple  of 

(35) ,  obtaining 

2>T„yst  *  b'*)  -7T-V  *  -(2't/y ’)  v*j: 
:(h'-x'v'V')/y'  .  (3?) 

In  orthogonal-function  space  this  is  a  set  of  seven  linear  algebraic 

equations  in  3  T  /  St  and  TC.  ,  ...,  %.  ,  whose  coefficients  are 

o  •  6 

obtained  by  transforming  X’/Y’,  T  ,  and  Z'T/Y’  to  orthogonal-function 
space,  and  whose  right-hand  sides  are  obtained  by  transforming 
(H1  -  X'  V  F')/Y'  .  Once  these  algebraic  equations  have  been  solved, 
the  integration  of  the  model  is  straightforward. 


5.  Sources  anu  sinks 

As  we  have  already  noted,  the  vertical  averages  of  F,  H,  and  (. 

represent  vertical  exchanges  of  momentum,  sensible  and  latent  heat,  and 
water.  We  shall  assume  that  the  exchanges  across  the  ocean-atmosphere 
interface  are  proportional  to  the  differences  of  appropriate  quantities 
across  this  interface,  with  the  same  factor  of  proportionality  k  .  We 
could  appeal  to  Ekman-layer  theory  to  determine  a  suitable  value  for  k, 
but,  in  view  of  the  drastic  s implf f feat  ions  already  introduced,  we  can 
hardly  justify  anything  more  involved  than  simply  choosing  a  time 

period,  say  five  days,  for  1/k  .  Exchange  of  momentum  across  the 

i 

surface  p  =  ^  ,  where  there  Is  no  Ekman  layer,  will  be  made 
proportional  to  the  shear,  with  a  damping  coefficient  k'  considerably 
smaller  than  k 

To  formulate  precipiatlon  we  shall  assume  that  clouds  have  a  "half 
life",  i.e.  ,  that  during  a  certain  time  interval  the  clouds  give  up  a 
certain  fraction  of  their  water  as  rain.  The  damping  time  1/Jj  should 

be  considerably  shorter  than  1/k  ,  perhaps  one  day.  Denoting  the 

radiative  contributions  to  H  and  E  by  Hj.  and  Ef.  ,  we  have 


f  =  -w  vV 


(38) 


F 


•  * 

-W  V  T 


(39) 


H/c^s  -Vi(t-3  +  Hjtr 
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6-  =  -  h(^-s')  -  ,  M) 

E/t:-ck(^TtlL/cr)(s-ri)t£r/c  ,  <«> 

where  c"  =  p  c  /tfp  -  p  )c).  The  factor  V/  X  appears  in  the 

O  ^  ®  ' 

precipitation  term  in  (41)  because  the  vertical  average  of  w  is  (^/X)w. 

Rariiation  is  traditionally  a  complicated  process,  and,  despite  our 
efforts  at  simplification,  our  formulation  will  reflect  this  fact.  Even 
if  we  are  interested  simply  in  an  atmosphere  with  an  absorbing 
constituent,  the  appropriate  radiation  formulas  will  be  highly  dependent 
upon  what  tills  constituent  is,  and  our  radiation  formulas,  to  a  greater 
extent  than  our  other  formulas,  will  be  based  upon  the  supposition  that 
it  is  water. 

We  shall  express  the  incoming  solar  radiation  in  terms  of  a 
planetary  temperature  T ^  ;  this  is  the  temperature  which  a  black  body 

covering  the  entire  sky  would  have  to  have,  in  order  to  be  equivalent  to 
tiie  sun.  We  shall  allow  some  of  the  incoming  solar  radiation  to  be 
reflected  by  clouds;  the  remainder  will  pass  through  the  atmosphere  and 
heat  the  ocean.  The  fractional  cloud  cover  a  ,  which  to  a  first 
approximation  will  equal  the  reflected  fraction  of  the  solar  radiation, 
vi 1 l  be  a  function  of  the  relative  humidity.  This  appears  preferable  to 
letting  the  cloud  cover  depend  upon  total  water  content;  we  are  not 
aware,  for  example,  that  cloudiness  is  greater  in  the  tropics  than  in 
the  polar  regions.  A  reasonably  suitable  function  is  a  =  (v/u)4. 


A 


The  ocean  will  in  turn  emit  lonj;-wave  radiation  as  a  black  bod'. 
The  cloud-free  portion  of  the  atmosphere  will  be  assumed  to  possess  a 
spectral  window,  through  which  a  fixed  fraction  a'  of  the  long-wave 
radiation  will  pass;  the  remaining  long-wave  radiation  will  be  partially 
absorbed  and  reemitted  by  atmospheric  water  vapor.  The  cloudy  portion 
of  the  atmosphere  will  behave  like  the  cloud-tree  portion,  except  that 
there  will  be  no  spectral  window. 

Our  expressions  for  emission  of  long-wave  radiation  will  need  to 

contain  the  atmosphere's  fractional  emisslvfty,  and  the  temperature  at 

which  the  emission  occurs.  In  a  model  with  such  low  vertical 

resolution,  we  balk  at  considering  the  radiative  transfer  from  level  to 

level,  and  find  no  reason  for  using  anything  more  complicated  than 

Simpson's  method.  Basically,  Sitnpson  (1928)  treated  the  atmosphere  as 

completely  transparent  In  the  8.5-11  micron  band  (the  window),  while 

below  7  and  above  14  microns  he  treated  any  layer  containing  0.3  kg  of 
o 

water  per  m  of  cross  section  as  completely  opaque.  Simpson's 
atmosphere  would  therefore  radiate  to  space  in  the  latter  wave  lengths 
at  the  temperature  of  the  uppermost  such  layer.  Between  7  and  8.5,  and 
between  11  and  14  microns,  more  water  vapor  would  be  needed  to  render 
the  atmosphere  opaque. 

Paraphrasing  Simpson's  treatment,  we  shall  as  a  first  approximation 
allow  the  cloud-free  portion  of  the  atmosphere  to  radiate  upward  and 
downward,  at  temperatures  T  '  and  T"  respectively,  with  the  fraction 
1  -  a'  of  the  intensity  of  black-body  radiation,  while  it  will  absorb 
the  fraction  1  -  a'  of  the  radiation  it  receives  from  the  ocean.  We 
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shall  neglect  any  variations  of  a'  which  ought  to  occur  with 
temperature.  We  shall  treat  the  cloud-covered  portion  of  the  atmosphere 
similarly,  except  that  a'  will  be  replaced  by  0  .  The  temperatures 

T'  and  T"  will  occur  at  pressures  p'  and  p"  ,  which  will  be  the 
pressure  levels  above  which  and  below  which  the  amount  of  water  vapor 
is  V*/2  ,  where  V*  =  0.3  kg/m2  .  Thus  p'  and  p"  are  supposed  to 
represent  the  centers  of  the  uppermost  and  lowermost  layers  whose  water 
vapor  content  is  V*  .  Letting  v*  be  the  value  which  v  would 
possess  if  the  water  vapor  content  of  an  entire  column  were  V*  ,  i.e., 


v*  -/>>  gV*/po  •  where  g  is  the  acceleration  of  gravity,  we  obtain, 
for  our  first  approximation. 


(43) 


(44) 


Kqs.  (41)  and  (44)  are  fairly  satisfactory  if  v/v*  is  large,  but 
if  v  <■  v  /2  ,  they  place  the  uppermost  layer  with  a  water  content  of 
V*  h?-l?w  Che  lowermost  layer,  while  in  actuality  such  layers  do  not 
exist  at  all.  We  shall  therefore  modify  (43)  and  (44)  by  replacing  v 
by  v  +  v*  .  Introducing  the  quantity  v'  =  v/(v  +  v* )  ,  which  is  near 


unity  when  v  is  large  but  is  small  when  v  is  small,  we  find  that 


(45) 


(46  ) 


r 
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whence 


I 

T  - 

I* 

T  - 


(4  7) 


( 4  K ) 


Eqs.  (45)  and  (46)  are  almost  Identical  to  (43)  and  (44)  when  v  is 
large,  but  they  make  p'  and  p"  approach  one  another  as  v  ->  U  . 

The  assumption  that  the  intensity  of  the  radiation  is  1  -  a' 

times  that  of  black-body  radiation  also  becomes  unrealistic  when  v  is 
small,  and  there  is  very  little  water  vapor  to  radiate.  Likewise,  a 
sufficiently  tenuous  cloud  layer  should  not  radiate  as  a  black  body.  We 
shall  adjust  for  this  situation  by  multiplying  the  emitted  and  absorbed 
radiation  by  v'  .  Finally,  a  tenuous  cloud  layer  shoud  not  be  a 
complete  reflector  of  solar  rdiaion,  and  we  shall  multiply  a  by  v' 
to  obtain  the  cloud  albedo. 

Collecting  our  results,  and  letting  a"  =  a  +  (1  -  a)a'  ,  we  find 

that 


ir/c  --  c"(^c;yr)  (- g"* t"*+  ) 


(49) 


(50) 


where  (f”  is  the  Stefan-Boltzman  constant 
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n.  Preliminary  numerical  runs 

We  have  developed  a  low-order  model  of  a  moist  circulation  for  the 
purpose  of  performing  some  comparative  numerical  experiments,  aimed  at 
clarifying  the  Importance  of  water  as  an  atmospheric  constituent.  We 
anticipate  comparing  pairs  of  runs;  the  members  of  a  single  pair  will  be 
alike  in  most  respects,  but  will  differ  in  the  value  of  one  or  possibly 
several  parameters.  Pairs  whose  members  differ  only  in  initial 
conditions,  and  which  may  therefore  be  expected  to  exhibit  similar 
statistical  properties,  can  be  used  as  controls.  Meanwhile,  as  is  often 
the  case  when  one  introduces  a  new  model,  our  first  studies  will  be 
directed  more  toward  toiling  us  something  about  the  model.  Inly  when  we 
feel  that  we  understand  the  capabilities  and  limitations  of  the  node] 
can  we  reasonably  expect  to  learn  something  about  the  atmosphere,  or 
about  atmosphere-tike  systems  which  the  model  attempts  to  simulate. 

A  necessary  prelude  to  performing  numerical  runs  is  the  selection 
of  numerical  values  for  the  various  constants.  We  have  chosen  our 
values  with  the  earth's  atmosphere  in  mind.  However,  to  reduce  the  time 

needed  for  transient  effects  associated  with  the  chosen  initial 

conditions  to  disappear,  we  have  made  the  oceanic  layer  interacting  with 
tlu  atmosphere  very  thin.  We  might  even  regard  the  underlying  surface 
as  wet  land  rather  than  ocean. 

The  dimensionless  values  of  f,  0,  po>  an^  ^  dre  unity.  We 

could  solve  the  equations  without  specifying  dimensional  values  at  all, 

but  in  that  case  it  would  be  difficult  to  compare  our  results  with  the 


real  atmosphere.  We  have  therefore  chosen  dimensional  values;  they  are 
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1/f  =  3  hours,  D  =  1830  km,  p0  •*  1000  mb,  and  R  =  287 

These  particular  values  have  the  advantage  of  making  the  unit  for 
temperature  exactly  100  degrees;  thus  no  further  conversions  beyond  a 
displacement  of  the  decimal  point  are  needed  for  atmospheric 
comparisons. 

Dimensionless  values  of  the  remaining  constants,  or,  where 
sufficient,  combinations  of  constants,  follow.  We  have  used  these 
values  in  our  numerical  runs  except  where  we  have  Indicated  otherwise. 


> 

0.  175 

y 

0.05 

X 

2/7 

L/cr 

25.0 

P 

0.3 

/ 

c 
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r 

0.  25 

k 

0.015 

k' 

0.001 

0.03 

c” 

1.0 

1 

CL 

0.5 

0.0006 

0.000061 
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In  our  first  runs  we  have  made  the  heat  received  from  the  sun,  us 
represented  by  T^,  Independent  of  geographical  location.  Our  interest 
in  these  runs  is  in  the  equilibrium  state  which  will  ultimately  be 
approached.  We  anticipate  that  in  this  state  T,  W,  and  S  will  be 
horizontally  uniform,  and  that  there  will  be  no  motion,  but  the  values 
of  T,  W,  and  S  which  should  correspond  to  a  particular  value  of  1, 
are  not  at  all  obvious.  The  problem  Is  therefore  non-trivial. 

In  the  final  Report  of  our  previous  contract  we  investigated  the 
same  problem.  We  discovered  to  our  surprise  that,  when  lay  in  the 

range  269  K-274  K,  there  were  two  stable  equilibrium  states,  and  also  an 
intermediate  unstable  equilibrium.  One  of  the  stable  equilibria  was 
warm,  and  one  was  cold,  relative  to  what  might  have  been  expected  with 
these  values  of  T^.  Outside  of  the  range,  high  values  of  produced 

only  warm  equilibria,  and  low  values  produced  only  cold  equilibria. 

Since  we  have  subsequently  made  a  number  of  changes  in  the  model,  it  is 
of  interest  to  see  whether  multiple  equilibria  still  occur. 

In  Fig.  1  the  heavy  curve  shows  the  equilibrium  values  ot  T  as  a 
function  of  Tq  .  We  see  that  there  are  still  multiple  equilibria. 

However,  the  range  of  values  of  T ^  for  which  such  equilibria  are  found 
has  been  narrowed  to  less  than  two  degrees.  One  might  say  that  the 
tendency  for  multiple  equilibria  to  occur  has  almost  disappeared. 

The  large  dots  on  the  curve  indicate  where  numerical  runs  were 
made.  The  dot  at  Tq  =*  275.  3,  T  =  277.  3,  representing  an  unstable 

equilibrium,  was  located  by  a  successive  approximation  procedure; 

previously  the  existence  of  the  unstable  equilibria  had  merely  been 
deduced  on  theoretical  grounds. 


Pig.  1.  Equilibrium  values  of  temperature  T  corresponding  to  assigned 
values  of  planetary  temperature  Tg  ,  as  given  by  the  model.  Heavy 
curve:  albedo  dependent  on  relative  humidity;  upper  thin  curve:  fixed 
low  albedo;  upper  heavy  curve  :  fixed  high  albedo. 
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In  the  earlier  report  we  hypothesized  that  the  existence  of 
multiple  equilibria,  and  the  accompanying  steep  slopes  of  the  curve  just 
outside  tne  multiple-equilibrium  range,  resulted  from  a  cloud-albedo 
feedback  process;  a  paucity  of  clouds  produces  a  low  albedo,  allowing 
more  solar  energy  to  be  absorbed,  thereby  evaporating  the  clouds  still 
more.  We  have  attempted  to  place  this  hypothesis  on  firmer  ground  by 
performing  some  auxiliary  runs  where  the  dependence  of  the  albedo  upon 

the  cloud  amount  has  been  eliminated.  Instead,  we  have  held  the  albedo 

4 

i  ixed,  in  one  set  of  runs  at  0.409b  =  0.8  ,  and  in  a  second  set  at 

4 

0.  ; 196  =  0.6  .  For  determining  the  long-wave  radiation,  the  dependence 

of  the  cloud  amount  upon  the  relative  humidity  has  been  retained. 

The  equilibrium  states  occurring  under  these  conditions  are 
indicated  by  the  two  thin  curves  in  Fig.  1,  the  upper  curve 
correspond! ng  to  the  lower  albedo.  The  tendency  for  multiple  equilibria 
hts  completely  disappeared,  and,  in  fact,  the  upper  curve  is  nearly  a 
str-tlght  line,  with  T  increasing  by  four  degrees  for  every 

three-degree  fncreso  in  T^  .  The  Importance  of  the  dependence  of  the 
albedo  upon  the  cloud  amount  is  thus  confirmed. 

Although  our  procedure  for  finding  equilibria  has  been  to  choose 
arbitrary  initial  conditions  and  let  the  model  run  until  equilibrium  is 
approached,  what  we  are  really  attempting  to  do  when  we  seek  ar. 
equilibrium  state  is  Co  solve  three  simultaneous  algebraic,  equations, 
two  of  which  are  highly  nonlinear.  These  equations  are  obtained  simply 
by  equating  li,  G,  and  E  to  zero.  The  simplest  of  these  equations, 

('.  =  0,  expresses  s  as  a  weighted  average  of  v  and  w;  with  the 


1 
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numerical  values  of  k  and  which  we  have  been  using,  s  *  (v+w)'2 

when  C  =  0.  As  we  noted  In  the  earlier  report,  t he  most  appropriate 

values  for  k  and  are  uncertain  in  any  case,  and,  by  decreasing  § 

by  a  factor  of  2,  we  reduce  the  equation  C  =  0  to  the  even  simpler 

equation  s  =  w,  which  implies  that  S  W.  This  result,  when  known, 
should  facilitate  the  analysis  and  interpretation  of  accompanying 
results.  Accordingly,  we  considered  it  desirable  to  examine  the 
equilibrium  curve  which  would  occur  when  ~  rather  than  0.03. 

The  curve  appears  as  the  right-hand  heavy  curve  in  Fig.  2.  The 

central  heavy  curve  is  the  original  equilibrium  curve,  copied  from  Fig. 
i.  The  range  of  values  of  T^  for  which  imi1tiple  equilibria  occur  has 
expanded  enormously,  and  the  difference  between  the  two  equilibrium 

values  of  T  ,  when  *  290K,  is  more  than  80  degrees.  One  suspects 

that  something  about  this  curve  is  unrealistic,  but,  in  any  case,  the 
multiple-equilibrium  phenomenon  is  seen  to  be  rather  sensitive  to  the 
chosen  numerical  valus  of  some  of  the  constants. 

The  contrast  between  the  curves  suggests  that  multiple  equilibria 

might  disappear  altogether  if  is  increased  rather  than  decreased. 

The  left-hand  heavy  curve  in  Fig.  2  has  been  constructed  from  runs  with 
=0.06.  There  are  indeed  no  multiple  equilibria;  there  are  only 
steep  slopes. 

The  variations  of  the  equilibrium  curve  are  in  qualitative 

agreement  with  what  one  might  expect  if  a  cloud-albedo  process  is 

operating.  The  low  value  of  ,  which  produces  the  widely  separated 


equilibrium  values,  corresponds  to  clouds  which  precipitate  their  water 


Fi\.  T. .  Foui  i  ii..j  i  um  values  of  temperature  T  corresponding  to  assigned  values 
if  nianetary  temperature  Tp  ,  as  given  by  the  model,  with  albedo  dependent 
01  relative  humidity,  for  various  choices  of  parameters  k  and  affectinq 
•vaporation  and  precipitation.  Left  neavy  curvet  medium  k  ,  high 
1  ”‘t  thin  curve:  low  k  ,  medium  ;  central  curvi' :  medium  k  , 

rv-uium  ;  right  thin  curve:  high  k  ,  medium  ;  right  heavy  curve: 

medium  k  ,  low 


only  with  some  difficulty,  and  which  therefore  tend  to  persist. 


The 


high  value  of  $  corresponds  to  clouds  which  are  rapidly  attenuated  by 
precipitation,  so  that  the  cloud-albedo  feedback  process  cannot  operate 
with  sufficient  strength  to  produce  the  multiple  equilibria. 

For  good  measure  we  have  made  two  more  sets  of  runs,  corresponding 
to  halved  and  doubled  values  of  k  ,  with  the  original  value  of  3 
The  results  appear  as  the  thin  curves  in  Fig.  2.  Qualitatively, 
decreasing  k  is  like  increasing  ^  ,  and  vice  versa,  hut  the  effect 
is  not  nearly  so  great.  Something  more  than  the  simple  ratio  3  /k  is 
involved;  the  heavy  and  the  thin  right-hand  curves,  which  are  far  from 
identical,  both  result  from  runs  with  3v/>  =  k,  and  so  with  S  =  W. 


Our 

final 

run 

allows  for  horizonal 

variat i ons 

and 

accompany  1  ng 

motions; 

it  has 

been 

made  with  Tq  =  277.0 

+  4.  0  $>  (in 

1 3 

°K). 

The  initial 

conditions  were 

obtained  by  superposing 

a  fairly  small 

perturbat ion 

At=  o,  AT=  Aw  =  AS  s  on  the  zonal  (but  not  steady)  state  'f'  =  n, 
T  =  W  =  S  =  T^,  Table  1  shows  the  development  of  the  temperature 
field,  expressed  in  terms  of  orthogonal  functions. 

Because  the  mean  value  of  T  ^  lies  close  to  the 

multiple-equilibrium  range,  the  low  and  high  latitudes,  where  T  ^ 
exceeds  or  falls  short  of  277.0  by  several  degrees,  correspond  to  very 
high  and  very  low  equilibrium  temperatures.  Consequently  T3  ,  which 
represents  the  south-north  temperature  contrast  and  the  accompanying 
westerly  wind  shear,  Increases  from  its  initial  value.  During  the  first 
month  the  waves,  represented  by  T, ,  T. ,  T^  ,  and  T^.  ,  damp  considerably 
and,  as  indicated  by  the  succession  of  positive  and  negative  values, 
progress  slowly  from  west  to  east. 


i, 
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Table  L.  Values  of  T  at  selected  times  In  a  particular 
integration  of  the  low-order  model,  with  =  277.0  +  4. 


t  (days) 

To 

T1 

T 

2 

T3 

T4 

T5 

0 

277.00 

1.00 

0.00 

4.00 

0.00 

0.00 

4 

276. 32 

0.  57 

0.56 

4.49 

0.00 

0.01 

8 

276.28 

-0.10 

0.65 

4.7  3 

-0.01 

0.00 

12 

276.32 

-0.46 

0.28 

4.93 

0.00 

-0.02 

16 

276.38 

-0.35 

-0.16 

5.16 

0.02 

0.00 

20 

276.44 

0.00 

-0.28 

5.38 

0.00 

0.02 

24 

276.51 

0.20 

-0.08 

5.58 

-0.02 

0.00 

28 

276.57 

0.11 

0.11 

5.80 

0.00 

-0.01 

32 

276.64 

-0.05 

0.09 

6.01 

0.00 

0.00 

36 

276.70 

-0.07 

-0.03 

6.22 

0.00 

0.00 

48 

2/6.90 

-0.01 

0.03 

6.86 

0.00 

0.00 

60 

277.09 

0.01 

0.01 

7.50 

0.00 

0.00 

72 

277.28 

0.02 

0.00 

8.15 

0.00 

0.00 

84 

27 7 . 46 

0.03 

0.04 

8.80 

c.oc 

0.00 

96 

277.60 

-0.25 

0.04 

9.43 

0.02 

0.00 

108 

2/7.63 

1  .70 

-0.78 

9.15 

0.04 

-0.14 

120 

2/2.6! 

1.40 

1.04 

7.34 

0.44 

0.18 

1  32 

270.98 

1.09 

0.22 

8.02 

0. 06 

0.00 

i  44 

270.07 

0.53 

0.88 

8.34 

0.04 

0.05 

1  56 

269.0) 

-0.68 

1.02 

8.41 

-0.C2 

0.15 

168 

267.66 

-1.33 

-0.  55 

8.37 

-0.16 

-0.05 

180 

266.  34 

0.01 

-1.23 

8.  -;6 

-0.05 

-0.12 

240 

260.86 

1 . 07 

-1.03 

8.43 

0.07 

-0.06 

300 

257.44 

1  .23 

0.59 

8 . 38 

-0.08 

-0.02 

160 

257.13 

0.  53 

0.71 

8.25 

-0.14 

-0.1 1 

420 

259.02 

-0.32 

0.36 

8.4'. 

0.04 

-0.11 

480 

260. 77 

0.37 

0.31 

8.92 

-0.13 

-0.04 

600 

263.49 

0.22 

-0.40 

10.21 

-0.  04 

0.  1  7 

7  2  0 

265.43 

0.27 

0.42 

11.34 

-0.  I  6 

-0.1) 

840 

266.43 

-0.06 

-0.  34 

11.59 

-0.26 

0.4  3 

960 

265.00 

0.13 

-0.69 

10.10 

'1.11 

0.02 

1080 

265.00 

0.  51 

0. 16 

10.10 

-0.07 

-0.  2  1 

numeric a 

$3 

V°K) 


0.00 

0.01 

0.01 

0.02 

0.03 

0.03 

0.03 

0.03 

0.03 

0.03 
0.02 
0.01 
0.00 
0.01 
0.00 
0.23 
0.  96 
0.49 
0.40 
0.  34 
0.  14 

0.  14 
0.20 
0.27 
1  .  22 

1 . 63 

1  .  94 

2.44 
2 . 7  0 

2.63 
2.26 

2  .  26 


During  the  third  month,  at  about  the  time  when  reaches  8.5 

degrees,  the  flow  becomes  baroelinically  unstahle,  and  the  waves  are 

suddenly  rejuvenated.  The  waves  derive  their  energy  from  the  zona] 
energy  (kinetic  or  available  potential),  so  that  by  the  end  of  the 

fourth  month  has  fallen  below  its  critical  value,  and  the  waves  have 

ceased  growing.  For  the  remainder  of  the  first  year  Tj  ,  and  the 
amplitude  of  the  waves,  remain  relatively  constant. 

A  feature  of  the  earlier  runs  which  is  apparent  also  In  this  run  is 
the  long  time  which  is  sometimes  needed  for  transient  effects  to  die 

out.  With  the  rejuvenation  of  the  waves,  the  mean  temperature  I'0 

suddenly  drops,  and  continues,  until  by  the  end  of  one  year  it  has 

fallen  a  full  20  degrees.  Presumably  the  waves  have  something  to  do 
with  the  decrease  of  T0  ,  but  the  effect  is  not  what  one  would  expect 
in  a  dry  model. 

At  about  tills  time  another  phenomenon  appears;  ,  which  has 

been  rather  small,  becomes  appreciably  positive.  A  positive  value  of  Tfc 
together  with  a  positive  value  of  Tj  ,  means  that  the  latitude  of 
maximum  westerly  shear  has  shifted  from  the  center  of  the  channel  toward 
the  south.  Evidently  the  current  is  less  unstahle  baroelinically  In  its 
new  position;  before  the  end  of  the  second  year  Ts  lias  Increased  by 
three  degrees,  with  no  accompanying  intensification  of  the  waves.  The 
mean  temperature  has  also  recovered  some  o.  its  loss,  which  supports  the 
suggestion  that  the  waves  may  have  played  a  role  in  the  loss. 

Few  variations  are  evident  In  the  latter  half  of  the  third  year 


other  than  the  progression  of  the  waves;  in  fact,  the  day-to-day  values, 


H't  presented  in  Table  2,  show  almost  no  variations.  Apparent iv  the 
transient  effects  have  finally  disappeared. 

it  is  surprising,  and  disturbing  from  the  point  of  view  <j : 
computational  economy,  chat  years  are  required  for  the  types  of  change 
that  one  might  have  expected  to  occur  in  months.  However,  aside  f  rnir. 
cne  slow  decay  times,  the  model  appears  to  behave  reasonably;  wuvt s 
develop  when  the  zonal  flow  is  sufficiently  strong,  and  they  progress  in 
a  reasonable  manner.  We  believe  that  the  model  is  ready  f.,r 

initial  "production  runs". 


/.  Ci including  remarks 


As  we  have  noted,  our  intention  is  to  use  our 
form  to  Investigate  the  ’  mport -nu  r  ui  w  iti  r 
constituent.  Meanwhile,  it  is  not  inappropriate 
further  changes.  Even  though  we  may  t  eel  u:a'  v 
time  to  tiie  development  process,  there  remains  tin 
model,  when  put  to  use,  will  prove  to  be  1  ;u  >.  i  i 
element . 


Model  in  its  nrts.uii 
as  an  atmospheric 
to  uiit  i  v  i  pa  t  e  still 
have  devoted  enouj'n 
possibility  that  i.a 
lg  in  some  essential 


Perhaps  the  horizontal  resolution  will  prove  to  be  toe  low. 
Fortunately  wo  have  constructed  the  mode!  so  that  we  can  inerts  the 
number  of  orthogonal  functions  and  the  accompanying  number  oi  grid 
points  without  making  any  other  changes.  Perhaps  we  shall  have  to 


formulate  the  clouds  (n  a  different  manner 


We  believe,  however,  that  the  most  serious  shortcoming,  from  i 

theoretical  point  of  view,  is  the  assumption  of  a  fixed  lapse  rate 
temperature  with  elevation.  Amplifying  waves  should  limit  tie  ::  uUKati 
growth  by  making  the  flow  on  which  tliey  are  superposed  less  unstable 
baroclinically.  One  way  in  which  they  can  do  this  is  by  reducing  the 
north-south  temperature  contrast;  our  model  allows  them  to  do  tills. 
Anther  way  is  by  transporting  heat  from  low  to  high  elevations,  thereby 
stabilizing  the  lapse  rate.  A  model  admitting  this  possibility  would  be 
more  realistic. 

Or.e  possible  scheme,  which  is  the  one  which  we  used  in  our 

vacillation  study  with  a  dry  model,  is  to  allow  the  lapse  rate  to  vary 
with  time  but  not  spatially.  To  make  our  present  model  capable  if  doing 
this,  i.e.,  to  make  X  a  function  of  time  alone  Instead  of  a  constant, 
we  would  need  to  make  some  assumption  concerning  the  transfer  of  energy 
from  one  level  in  the  atmosphere  to  anotiier  by  radiation  and  other 
processes  -  an  assumption  which  we  have  been  expressly  trying  to  avoid. 
Also,  in  our  present  model  we  have  stored  a  number  of  functions  of  N- 
as  constants;  these  would  become  variables.  Nevertheless,  we  see  no 

insurmountable  di f f Icul t ies. 

As  already  stated,  we  believe  that  our  immediate  task  Is  to  make 
some  production  runs.  Eventually,  however,  we  believe  that  we  should 
repeat  some  of  these  runs  with  a  model  with  a  variable  X  ,  to  see 


whether  the  conclusions  continue  to  stand 
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